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Abstract We present broad-band imaging with the Subaru Telescope of a 25'x25' field surrounding the radio galaxy TN J1338~ 
1942 at redshift z = 4. 1 . The field contains excesses of Lyman-a emitters (LAEs) and Lyman break galaxies (LBGs) identified 
with a protocluster surrounding the radio galaxy. Our new wide-field images provide information about the boundary of the 
protocluster and its surroundings. There are 874 candidate LBGs within our field, having redshifts in the range z = 3.5-4.5. An 
examination of the brightest of these (with i' < 25.0) shows that the most prominent concentration coincides with the previously 
discovered protocluster The diameter of this galaxy overdensity corresponds to ~ 2 Mpc at z = 4, consistent with the previous 
estimation using LAEs. Several other concentrations of LBGs are observed in the field, some of which may well be physically 
connected with the z = 4.1 protocluster The observed structure in the smoothed LBG distribution can be explained as the 
projection of large-scale structure, within the redshift range z = 3.5 - 4.5, comprising compact overdensities and prominent 
larger voids. If the 5-8 observed compact overdensities are associated with protoclusters, the observed protocluster volume 
density is ~ 5 x 10"* Mpc"^, similar to the volume density of rich clusters in the local Universe. 

Key words. cosmology:observations — early universe — galaxies:clusters:general — galaxies:high-redshift — large-scale 
structure of universe 

1 . Introduction ters may be typical. Because radio-sources are relatively short- 

lived (~ 10^ years; Blundell & Rawlings ll999> . the statistics 



^ . There is considerable evidence for galaxy overdensities athigh consistent with the progenitor of every rich local cluster 

■ ■ ' redshifts (z > 2; e^teidel et al. 1998, M0ller & Fynbo^mi having harboured a luminous radio galaxy at some stage in its 

ShimasakuetalEnnH Palunasetal. 2004, Ouchietal. 2005). existence 
In many cases these overdensities have been presumed to be 

associated with the ancestors of rich local clusters. At z > 2 Using the VLT, Venemans et al. C002^ spectroscopically 

the Universe is < 3 Gyr old, too short for these structures to confirmed 20 Lyman-a emitters (LAEs) in a 7' x7' field around 

have virialized (e.g. Venemans 2005 ). Hence these structures "zRG TN J1338-1942 at a redshift of z = 4.1 and identified 

are often called protoclusters. Most searches for protoclusters *ese LAEs with a z = 4.1 protocluster. Further evidence that 

have been limited by relatively small fields (typically smaller *is LAE overdensity was indeed associated with a protocluster 

than 10' X 10') provided by observations with the HST/ACS that revealed 

. „ . raj- 1 • 1 • 1 J an excess and non-uniform distribution of candidate Lyman 

An efficient way of finding protoclusters IS to use high red- , , , . , „^ ^ , , ,. , , 

, .„ ,. , . . , break galaxies (LBGs) around the radio galaxy (Miley et al. 

shift radio galaxies (HzRGs; e.g. Rottgering et al. I1994t as rr„„ . . , , 

tracers (Venemans et al. |2002l Kurk et al. Iffl. HzRGs are ^^^^ Overzier et al. 2006 ). 

large massive objects with many of the properties expected of Although the LAE search around TN J133 8-194 2 was ex- 
forming dominant cluster (cD) galaxies (West'T994). Although tended with a second 7' x 7' field (Venemans 2005 ), this was 
most protoclusters are not radio-loud, radio-selected protoclus- insufficient to determine the boundary of the protocluster. Here 

we present the results of a multi-color study of candidate LBGs 

* Based on data collected at Subai'u Telescope, which is operated from a 25' x25' region surrounding TN J1338-1942. The large 

by the National Astronomical Observatory of Japan. field-of-view (FOV) facilitates searches for LBGs out to the 
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boundary of the protocluster structure and beyond. The data 
also provide new information about large-scale structure and 
voids at z ~ 4. 

Throughout this paper, we use AB-magnitudes, Icr errors 
and adopt a flat, A-dominated cosmology with Q.^ = 0.3, Q.a = 
0.7 and Ho = 100 /i km s"' Mpc"' with h = 0.7. 

2. Data reduction and sample selection 

Deep multi-color imaging of the TN J1338-1942 field was car- 
ried out using the Subaru/Suprime-Cam instrument (Miyazaki 
I2002t on January 3P' and February P', 2003. Data reduction 
on the B-, Rq- and /'-band images was performed in a man- 
ner similar to that of the SDF and SXDF fields (Ouchi et al. 
I200ni2004ail . The FOV was 24.7' x 24.2' and the seeing had an 
equivalent FWHM of 0.98". Source extraction and photometry 
was done using SExtractor (Bertin & Arnouts 1996). Fluxes 
were measured in a 2" circular aperture with 3cr limiting mag- 
nitudes of 27.2, 26.9 and 26.5, respectively. Extending the ob- 
served power-law estimation for bright source counts (e.g. see 
Palunas et al. 2004), the initial object sample was limited to 
i' < 26.5 to obtain a photometric completeness of ~ 72% for 
the highest magnitude bin (A/' = 0.5). 

LBGs in an approximate redshift range z = 3.5 - 4.5 were 
selected using color selection criteria by Ouchi et al. J2004at . 
resulting in a sample of 874 LBGs. Based on the same work 
by Ouchi et al., the estimated contamination by interlopers and 
stars is ~ 6%, while the estimated completeness distribution 
function over redshift has an approximate gaussian function 
shape with a FWHM of 0.8, centered at z = 4 with a peak value 
of ~ 45%. In addition, a bright sub-sample of LBGs was con- 
structed having /' < 25.0. For this sub-sample of 125 objects, 
the galaxy colors at the corresponding LBG redshift range are 
well constrained. The completeness distribution function for 
this sub-sample has a FWHM of 0.8, centered at z = 4.1 with 
a peak value of ~ 90%, while the estimated contamination is 
< 1%. 

3. Analysis 

3.1. Projected density distribution of bright LBGs 

We first investigated the projected distribution of the bright- 
est LBGs by smoothing the spatial distribution of our bright 
sub-sample with a gaussian kernel. Structure identification is 
dependent on the size of the smoothing kernel. The FWHM 
of 5' was chosen to match the average distance between these 
LBGs, thereby optimizing the contrast between overdense and 
underdense regions. This FWHM is also similar to the angular 
size of LAE proto-clusters at z = 3 - 6 (Shimasaku et al. 120031 
Ouchi et al. 2005 Venemans et al. 2005 1, which improves the 
chance of detecting such structures in the LBG distribution. 
The smoothed LBG map was divided by a normalisation map 
to correct for undetected LBGs behind foreground objects and 
FOV boundaries, which causes some dense areas near the bor- 
der to be overemphasized. The result is shown in Fig.Q] 

The most significant projected overdensity in this bright 
LBG structure map (away from the border) contains the ra- 
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Figure 1. Projected distribution of 125 bright LBGs (open cir- 
cles; diameter scales with brightness; 21.5 < /' < 25.0) in the 
TN J1338-1942 field, including TN J1338-1942 (filled circle 
within star). The projected density contours (curved lines), ob- 
tained by gaussian smoothing, reveal overdense (solid lines; 
A - 0.25, 0.50, 0.75, 1.00 from edge to center) and under- 
dense regions (dotted lines; A = -0.25, -0.50, -0.75 from 
edge to center) relative to a mean density of 0.21 per square 
arcminute (dashed line). TN J1338-1942 inhabits a signifi- 
cant overdense area, probably associated with the protoclus- 
ter found by Venemans et al. (2002i. The rectangles represent 
the two fields that were used by Venemans et al. (I2002> and 
Venemans ( 2005 ) to search for LAEs. 



dio galaxy TN J1338-1942. We associate this overdensity 
with the protocluster at z - 4.1 previously discovered in the 
smaller fields accessed by the VLT (LAEs in two 7' x 7' fields; 
Venemans et al. 2002, Venemans 2005 1 and ACS/HST (LBGs 
in a 3' X 3' field; Miley et al. (20041 Overzier et al. 2006 ). Taking 
A = 0.5 as the boundary (with A - (L - (S))/(Z), where Z and 
(E) are the local and average projected LBG densities, respec- 
tively), (i) the diameter of the overdensity is ~ 5', correspond- 
ing to a (proper) size of ~ 2 Mpc at z = 4 and (ii) the location 
of the radio galaxy is in the western part of the structure. The 
protocluster size and the relative location of the radio galaxy 
are similar to those found by Venemans et al. for LAEs. 

3.2. LBG overdensity in redshift space 

We determined the overlap between the whole LBG sample 
and the LAE sample by Venemans et al. (I2002> and Venemans 
(2005 ) to obtain spectroscopic redshifts for several LBGs. 
There were 86 LAE candidates found in a narrow redshift range 
Az = 4.087 - 4.119 surrounding TN J1338-1942 (fields plot- 
ted in Figs, n] and Spectroscopic confirmation of redshift 
followed for 38 LAEs, including TN J1338-1942. 
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Within the same area, we identified 104 candidate LBGs 
(out of the whole LBG sample). Of these LBGs, 7 are also 
spectroscopically confirmed LAEs, thus obtaining a redshift for 
these 7 LBGs. Based on the contamination fraction, we expect 
6 of the 104 candidate LBGs to be interlopers. To see whether 
the 7 confirmed LBGs represent a significant overdensity, we 
estimated the expected number of LBGs in the redshift range 
Az using Monte-Carlo simulations. We randomly assigned red- 
shifts to 98 objects, using the completeness distribution func- 
tion as redshift distribution, and counted the number of objects 
in Az. We repeated this procedure 10,000 times and found that 
the expected number of LBGs in redshift range Az is 3 .6 + 1 .9. 
As a result, the 7 spectroscopically confirmed LBGs might in- 
dicate a modest LBG volume overdensity of 6 - 1.0+1.0 
(with (5 = (p - (p))/(p), where p and (p) are the local and av- 
erage LBG volume density, respectively) in close vicinity of 
TNJ1338-1942. 

The estimated overdensity above is a lower limit, because 
not all high redshift galaxies have Ly-o- emission. Steidel et al. 
(20001 found that at z = 3 only 20 - 25% of the LBGs satisfy 
typical LAE selection criteria. Assuming that this is also true at 
Z = 4, this implies that the 7 confirmed LBGs represent a true 
number of at least 28 LBGs within Az. As these are expected to 
be part of the 104 observed LBGs, the LBG volume overdensity 
within Az is increased to at least 6-1+4. Note that still 
~ 70% of the LBGs lie outside Az. This excess of LBGs in 
redshift space is consistent with the location of the protocluster 
that harbours TN J1338-1942. 

3.3. Projected density distribution of all LBGs 

Using the same technique as in Sec. 13.11 (but using a FWHM 
of 2' to match the average distance between the LBGs), the 
whole sample of 874 LBGs was used to make a second struc- 
ture map, which is shown in Fig.|5] In this map, TN J1338- 
1942 also inhabits a clear but less prominent projected overden- 
sity of galaxies. In addition to this overdensity associated with 
the previously known protocluster, 4-7 other intriguing peaks 
are seen in the large-scale structure distribution, surrounded by 
larger regions of relatively empty space. These overdensities 
may be associated with protoclusters within the redshift range 
z = 3.5 - 4.5, while the underdense regions indicate the pres- 
ence of large voids. The typical transverse size of the overden- 
sities is ~ 2 Mpc at z = 4, while the underdensities are more 
than twice this size. Without spectroscopic data, we cannot es- 
tablish whether the other overdensities are physically hnked to 
the protocluster. 

The overdensity surrounding TN J 1338-1 942 contains rel- 
atively many bright LBGs compared to the other overdensi- 
ties, which suggests that the TN J1338-1942 protocluster is the 
most massive structure within the observed volume. The strong 
clustering of bright LBGs at one particular position within the 
FOV agrees with the observation that at z ~ 4, brighter LBGs 
have larger clustering lengths than fainter LBGs (Allen et al. 
l200?l Ouchi et aI. I2()()4bt . 

Similar structure maps were created for 20 mock samples of 
874 random points each with the same positional constraints as 
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Figure!. Projected distribution of 874 LBGs with 21.5 </' < 
26.5 in the TN J1338-1942 field, similar to Fig.C] The con- 
tours are relative to a mean density of 1 .47 per square arcminute 
(dashed Une). Like in Fig.^ TN J 1338-1 942 inhabits a signif- 
icant overdense area, probably associated with the protocluster 
found by Venemans et al. (i2002j . 



the LBG sample. Visual comparison between the detailed LBG 
map and the mock maps shows that for the latter, the overdensi- 
ties are larger but lower in amplitude, while the underdensities 
are smaller and more isolated. Basically, these over- and under- 
densities have similar sizes and amplitudes, very different from 
what is observed in the LBG map. 

Figure |3] shows the projected density distribution function 
(PDDF) of the detailed LBG map and the mean PDDF of the 
mock maps. Using a Kolmogorov-Smirnov test (e.g. Press et 
al. '1992^, we found that the probability that the PDDF of the 
LBG map is drawn from an underlying distribution equal to 
the mean PDDF of the mock maps is 43%. For the 20 mock 
maps, the probability that they are drawn from the mean PDDF 
is much higher (> 99.9%). In Figure |3 it can be seen that for 
the LBG sample there is significantly (> 3cr) more area with 
-0.75 < A < -0.5 and A > 0.65 than for the mock samples. 
This is consistent with the underdensities (presumably voids) 
being larger than the more strongly peaked overdensities (pre- 
sumably protoclusters). 

3.4. LBG angular and spatial correlation 

The two-point angular correlation function (ACE; Peebles 
119731 11980I I of the whole LBG sample was calculated, using 
estimators by Landy & Szalay (1993) and Hamilton ( I1993> . 
which gave practically identical results. The estimators are neg- 
atively off^set from the true ACE due to the diff'erence between 
the measured average and the true average LBG density in the 
restricted FOV (the 'integral constraint'). When assuming a 
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Figures. Plotted here are the (normalised) projected density 
distribution function (PDDF) of the whole sample of 874 LBGs 
(solid line) and the combined PDDF of 20 equally sized mock 
samples (dashed line is mean, dotted lines are mean + 3cr, 
where <t is the standard deviation between the mock samples). 
At several intervals, the PDDF of the LBG sample deviates 
from the mock samples by more than 3cr (thick solid line). The 
smoothed LBG distribution has relatively less area with a den- 
sity close to mean, relatively more area out to higher densities 
and relatively more area with lower densities, which character- 
izes the presence of compact overdensities with high peaks and 
extended underdensities. 



power-law form a)(ff) = 0'^ for the true ACF, the offset can 
be estimated following Roche et al. ( 2002). Iteratively fitting a 
power-law and estimating the offset converged to solutions for 
both a variable {J3 - \.\ and a fixed slope (fi - 0.8; e.g. 

Peebles "^1980 1 power-law. We found that for all cases the clus- 
tering amplitude is significantly larger than its uncertainty 
(> 6cr), confirming that there is a non-random clustering signal 
present in our LBG distribution. 

After correcting the clustering amplitude for contamina- 
tion (Ouchi et al. 2004bl. the inverse Limber transformation 
(Efstathiou et al. 199ll was used to calculate the (comov- 
ing) spatial correlation length. We found correlation lengths 
of ro = (3.7 + 0.7) h'^ Mpc and ro = (4.6 + 0.4) /i"' Mpc 
for the variable and fixed slope fit, respectively. These results 
are similar to those found by Ouchi et al. ( 2004b ), indicating 
that (within the observed volume) the clustering properties of 
LBGs in the radio galaxy field are not significantly different 
from blank fields. 

3.5. LBG void probability function 

A disadvantage of using the two-point ACF for detecting clus- 
tering in galaxy distributions is that it is less sensitive to non- 
gaussian density fluctuations. One can use n-point ACFs to im- 
prove the detection of the latter, but error-bars becomes exceed- 
ingly large with higher n when the number of galaxies is fixed. 
A complementary approach is to focus on voids instead. The 
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Figure 4. The two-dimensional void probability function 
(VPF) of the LBG sample (filled dots; error bar is poissonian 
error), compared with the mean VPF of 20 mock samples (open 
squares; error bar is combination of poissonian error and stan- 
dard deviation between mock catalogs) and the theoretical VPF 
for poissonian distributions (solid line). For radii 6 > 40", the 
underdensities in the projected LBG distribution are signifi- 
cantly larger than the underdensities in mock and poissonian 
distributions. 



two-dimensional void probability function (VPF; White [T9791 
defines the fraction of circular areas at random positions in the 
FOV which contain no galaxies. The VPF was calculated for 
both the whole LBG sample and the 20 mock samples from 
Sec. l3.3l and plotted in Fig.|4] The mean VPF of the mock sam- 
ples is very similar to the theoretical VPF for poissonian dis- 
tributions (which is P\/pp(6) = exp [-.7r0^(2)]). At larger radii, 
the VPF of the full LBG sample is significantly higher than 
both the mean mock VPF and the theoretical VPF, meaning 
that the underdensities in the projected LBG distribution are 
relatively large compared to underdensities found in random 
distributions. This result is similar to the result of Palunas et al. 
(^2004 1 at z = 2.34, using a sample of 34 LAEs in a FOV similar 
to ours. 

3.6. Protocluster volume density 

The comoving volume that is observed using the LBG selec- 
tion criteria by Ouchi et al. ( 2004a i has a size of approxi- 
mately L8 X lO*" Mpc-'. Assuming that 5-8 observed over- 
densities in the LBG map are indeed associated with protoclus- 
ters, the volume density of protoclusters within z = 3.5 - 4.5 is 
-5x10*' Mpc ^ This agrees with the estimate of Venemans 
J2005I . who found a density of ~ 6 x 10 Mpc"^ for LAE 
protoclusters at z - 2-5.2. They also report on similar re- 
sults from Steidel et al. (fT998t . based on LBGs at z = 2.7 - 3.4 
(3 X 10"* Mpc"^), Shimasaku et al. (2003 1, based on LAEs at 
z = 4.9 (7 X 10"'' Mpc"^^) and Ouchi et al. (2005,1, based on 
LAEs at z = 5.7 (> 2 x 10"'' Mpc"^^). 

Our estimated volume density of protoclusters can be com- 
pared with the volume density of local rich clusters. Values 
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found for rich cluster density at low redshift (z < 0.1) lie 
in the range (2 - 4) x 10^'' Mpc^^ (Bahcall & Soneira [MSl 
Postman et al. T992 Peacock & West rT992l Zabludoff et al. 
1199 3 Mazure et al. ■1996j . Our result is consistent with this 
number range, providing further evidence that the observed 
overdensities at high redshift are progenitors of rich clusters 
in the local Universe. 

4. Conclusions 

We draw several conclusions from the present observations: 

(i) TN J1338-1942 is located in an overdensity of LBGs, both 
in projection and in redshift space. The new wide-field results 
are consistent with previous observations which revealed the 
presence of a protocluster through the overdensities of LAEs 
(Venemans et al. 2002i Venemans l2005t and LBGs (Miley et 
al. p004 Overzier et al. I2006t . This further supports the hy- 
pothesis that HzRGs are located in dense environments. The 
apparent size of the overdensity and the relative position of 
TN J1338-1942 within the overdensity are similar to that found 
by Venemans et al. for LAEs. 

(ii) There are 4 - 7 additional overdensities in the projected 
LBG distibution, similar to the one harbouring TN J 1338- 
1942. These may well be due to protoclusters at z = 3.5 - 4.5 
and one or more of these overdensities could well be physically 
related to the TN J 1338-1 942 protocluster. 

(iii) The spatial distribution of our complete LBG sample is 
consistent with a Universe at z ~ 4 that comprises a web 
of compact galaxy overdensities (protoclusters) embedded in 
larger regions of galaxy underdensities (voids). The statistics 
of the overdensities are consistent with the local volume den- 
sity of rich clusters. 

Spectroscopic measurements are needed to investigate 
whether there is a physical connection between some of the 
outlying observed galaxy overdensities and the overdensity 
corresponding to the TN J1338-1942 protocluster Such obser- 
vations could enable the cosmic web at z = 4. 1 to be traced over 
distances of tens of megaparsec. Furthermore, similar measure- 
ments on other z > 2 protoclusters would be useful for con- 
straining the development of large-scale structure in the early 
Universe. 
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